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T— I Abstract We consider networks with two types of nodes. The w-nodes, 

called centers, are hyperconnected and interact one to another via many u- 

y-^ nodes, called satellites. This centralized architecture, widespread in gene net- 

rH works, possesses two fundamental properties. Namely, this organization creates 

^ ^^ feedback loops that are capable to generate practically any prescribed pattern- 

(«-! ing dynamics, chaotic or periodic, or having a number of equilibrium states. 

• ^H Moreover, this organization is robust with respect to random perturbations of 

^~ the system. 

1 Introduction 

JI^ Flexibility and robustness are important properties of biological systems. Flex- 

QQ ibility means the capacity to adapt with respect to changes of environment 

T-H whereas robustness is the capacity to support homeostasis in spite of environ- 

^^ mental changes. Intriguingly, it seems that biological systems could be in the 

^^ same time robust and flexible. Development of an organism is robust to vari- 

^^ ations of initial conditions and environment, species can diversify in order to 

Oh better satisfy constraints imposed by a varying environment. 

. . We discuss here flexibility and robustness problems for gcnetical networks 

^ of a special topological structure as a model for flexible and robust systems. 

In these networks, highly connected hubs play the role of organizing centers 

(centralized networks). The hubs receive and dispatch interactions. Each center 

C^ interacts with many weakly connected nodes (satellites). Similar ideas, that 

such a "bow-tie" connectivity can play a role in robustness, have been proposed 
by (Zhao et al. 2006). In the field of random boolean networks (Kauffman 
1969), the phase transitions from chaotic to frozen (robust) phases were related 
to scale- freeness and heterogeneity of the network by (Aldana 2003). We show 
that centralized network are capable to produce a number of patterns, while 
being protected against environment fiuctuations. 

Network models usually involve interactions between transcription factors 
(TFs) (Rcinitz et al. 1991). In the last years, a great attention has been 
focused on microRNAs (He and Hannon 2004, Bartel 2009, Hendrikscon et al. 
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2009, Ihui et al. 2010). MicroRNAs (miRNAs) are short ribonucleic acid (RNA) 
molecules, on average only 22 nucleotides long and are found in all eukaryotic 
cells. miRNAs are post-transcriptional regulators that bind to complementary 
sequences on target messenger RNA transcripts (mRNAs) and repress transla- 
tion or trigger mRNA cleavage and degradation. Thus, miRNAs have an impact 
on gene expression and it was shown recently that they contribute to canaliza- 
tion of development (Li et al. 2009). (Shalgi et al. 2007) shows the existence 
of many genes submitted to extensive miRNA regulation with many TF among 
these "target hubs". Without excluding other applications, we consider regu- 
lation of TF by miRNAs as a possible example of a centralized network. For 
this particular situation we generalize the TF network models (Reinitz et al. 
1991) to take into account miRNA satellites and the centralized architecture. 
The interaction between network nodes is defined by sigmoidal functions that 
can be defined by two parameters: the maximum rates of production r^, and 
sharpness constants Ki. Other important parameters, play a key role, namely, 
the degradation constants A^ > of centers and satellites. 

We obtain a fundamental relation between the main network parameters. 
This relation ensures maximal robustness of the network with respect to random 
internal and external fluctuations, given a certain amount of flexibility defined 
as the number of attractors that are accessible to the network dynamics. Our 
mathematical results have a transparent biological interpretation: centralized 
motifs can be simultaneously flexible and robust. One can expect that miRNA 
molecules, being smaller with respect to TF, are more mobile and react faster 
to perturbations. This property plays a key role in the flexible and robust 
functioning of centralized motifs. 

The paper is organized as follows. Centralized networks are introduced in 
Section 2. We also formulate here an important assertion on the flexibility of 
general centralized networks. We show that these networks are capable to gener- 
ate practically all dynamics, chaotic or periodic, with any number of equilibrium 
states. To study robustness with respect to random fluctuations, in Section 3 
we consider a toy model of simple centralized TF - miRNA networks with a 
single center. We show here, by an elementary way, that centralized networks 
with mutually repressive hub-satellite interaction can produce many different 
robust patterns. 

2 Centralized networks 

Centralized networks have been empirically identified in molecular biology, where 
the centers can be, for example, transcription factors, while the satellite regu- 
lators can be small regulatory molecules such as microRNAs (Li et al. 2010). 
Notice that, in the last decades, the theory of so-called scale-free networks has 
become very popular. Scale- free networks (Barabasi and Albert 2002, Lesne 
2006) occur in many areas, in economics, biology and sociology. In the scale- 
free networks the probability P{k) that a node is connected with k neighbors, 
has the asymptotics Ck~'^, with 7 e (2,3). Such networks typically contain a 



few strongly connected nodes and a number of satellite nodes. Hence, scale-free 
networks are, in a sense, centralized. 

In order to model dynamics of centralized networks we adapt a gene circuit 
model proposed to describe early stages of Drosophila (fruit-fly) morphogenesis 
(Mjolness et al. 1991, Reinitz and Sharp 1995). To take into account the 
two types of the nodes, we use distinct variables Uj, Ui for the centers and the 
satellites. The real matrix entry Aij defines the intensity of the action of a center 
node j on a satellite node i. This action can be either a repression Aij < 
or an activation Aij > 0. Similarly, the matrices B and C define the action 
of the centers on the satellites and the satellites on the centers, respectively. 
Let us assume that a satellite can not act directly on another satellite (the 
principle of divide et impera). We also assume that satellites respond more 
rapidly to perturbations and are more diffusive/mobile than the centers. Both 
these assumptions are natural if we identify satellites as microRNAs. 

Let M, N be positive integers, and let A, B and C be matrices of the sizes 
N X M,M X M and M x N respectively. We denote by A^, Bj and Cj the rows 
of these matrices. To simplify formulas, we use the notation 

M M N 

y^ AijVj =- Ajt;, ^ Bjivi = BjV, ^ CjkUk = CjU. 
j=i 1=1 fe=i 

Then, the network model reads (we exclude diffusion effects): 

—j-^ficriAiV-hij-XiUi, (1) 

-^ = TjO (Bjf ^ C^u - hj) - XjVj. (2) 

We assume that the rate coefficients rj,ri are non- negative: ri,ri > 0. Here 
i = 1,...,N, j = 1,...,M and cr is a monotone and smooth (at least twice 
differentiable) sigmoidal function such that 

cr(-oo)=0, cr(+oo) = 1. (3) 

Typical examples can be given by the Fermi and Hill functions: 

'^(")=l + exp(-x)' '^^^"^^i^lT^' ^^) 

where Ka, p > are parameters and in the second case a; > 0. For x < we 
set anix) = 0. Analytical and computer simulation results are similar for both 
variants a and an- 

The parameters Xi,Xi are degradation coefficients, and hi, hi are thresholds 
for activation. 

Let us prove that the gene network dynamics defines a dissipative dynamics. 
In fact, there exists an absorbing set B defined by 

B = {w^{u,v) ■.()< vj < TjXj^, 0<u,< fjA~\ j = 1, ..., M, i = 1, ...,N}. 



One can show, by comparison principles for ordinary differential equations, that 

< Ui{x,t) < (j),{x) exp{~X,t) +t';XY^{1 -exp(-A,t)), 
< Vi{x, t) < (j)i{x) exp(-Aii) + r,X^ (1 - exp(-Ait)). 

Therefore, solutions of (II]), (pi) exist for all times t and they enters for the set 
S at a time moment Iq and then stays in this set for all t > Iq. So, our system 
defines a dissipative dynamics and all concentrations are positive if they are 
positive at the initial moment. In mathematical terms, the Cauchy problem 
(initial value problem) for our system is well posed. 

3 Complex dynamics of centralized networks 

Let us show that the centralized networks have a formidable power in dynamics 
generation. First, we will find an asymptotic simplification of the dynamics, 
then show that any dynamics, periodic, chaotic, or with a number of stable 
steady states can be approximated by centralized networks. 

3.1 Simplified dynamics when satellites are fast 

We suppose here that the u-variables are fast and the w-ones are slow. Then 
the fast u variables are slaved, for large times, by the slow v modes: one has 
u = U{v) + u, where m is a small correction. This means that, for large times, 
the satellite dynamics is defined almost completely by the center dynamics. 

To realize this approach, let us assume that the parameters of the system 
satisfy the following conditions: 

\A,l\,\Bul\C^J\,\h^l\h,\<Co, (6) 

where i = 1,2,. ..,N, i,l = 1,...,M, j = l,...,N, 

0<Ci<Aj, (7) 



and 

where 



n = uRi, f, = KRi, (8) 

\R,\,\R,\<C5, A, = kA,;, |A|<C6, (9) 

where k is a small parameter, and where all positive constants Ck are indepen- 
dent of K. 

Assertion 2.1. Under assumptions Im, [^, (M for sufficiently small k < 
AvQ solutions {u, v) of pp, M) satisfy 

u = U{v{t))+u{t), (10) 

where the j-th component Uj of U is defined by 

-XjUj = kGj{v), (11) 



where 

The function u satisfies estimates 

\u\ < cK^ + i?exp(-/3i), /3 > 0. (12) 

The V dynamics for large times t > Ci\ logKJ takes the form 

civ ' 

—^^KFi{u,v)+Wi, (13) 

where Wi satisfy 

\Wi\ < CK^ 

and 

F,{u, v) = R,a (B,v + C,C/(w) - h,) - X^v^. 

This assertion, known in computational biology as the quasi-steady state 
assumption, can be proved by well known methods from the theory of differential 
equations (Henry 1981). 

3.2 Realization of prescribed dynamics by networks 



Our next goal is to show that dynamics ( |13[ ) can realize, in a sense, arbitrary 
dynamics of the centers. To precise this, let us describe the method of realiza- 
tion of the vector fields for dissipative systems (proposed by Polacik 1991, for 
applications see, for example, Dancer - Polacik 1999, Rybakowski 1994, Vaku- 
lenko 2000). This method is based on the well developed theory of invariant 
and inertial manifolds, see Marion 1989, Mane 1977, Constantin, Foias, Nico- 
laenko and Temam, 1989, Chow-Lu 1988, Babin-Vishik 1988). One can show 
that there are systems enjoying the following properties: 

A These systems generate global semiBows S^ in an ambient phase space 
H . These semiBows depend on some parameters V (which could be elements of 
another parameter space B). They have global attractors and finite dimensional 
local attracting invariant C^ (continuously differentiable) - manifolds M., at 
least for some V. 

B Dynamics of Sip reduced on these invariant manifolds is, in a sense, "al- 
most completely controllable" . It can be described as follows. Assume the 
differential equations 

%=F{p), FeC\B-) (14) 

define a dynamical system in the unit ball i?" C R". 



For any prescribed dynamics {14) and any S > 0, we can choose suitable 
parameters V — V{n, F, 5) such that 

Bl The semiHow S\, has a C^- smooth locally attracting invariant manifold 
Ai-p diffeomorphic to the ball B"; 



B2 The reduced dynamics S^Im-p ^^ defined by equations 

= F{p,V), FeC\B^) (15) 



dt 



where the estimate 



\F - F\c^B^.) < 5 (16) 

holds. In other words, one can say that, byV, the dynamics can be speciGed to 
within an arbitrarily small error. 

Thus, all dynamics can occur as inertial forms of these systems. Such systems 
can be named maximally dynam.ically flexible, or, for brevity, MDF systems. 

Such dynamics can be chaotic. There is a rather wide broad in different def- 
initions of "chaos". In principle, one can use here any concept of chaos. If this 
chaos is stable under small C^ -perturbations this kind of chaos occurs in the dy- 
namics of MDF systems. To fix ideas, we use here, following Ruelle and Takens 
1971, Newhouse, Ruelle and Takens 1971 Smale 1980, Anosov 1995), such a def- 
inition. We say that a finite dimensional dynamics is chaotic if this generates a 
non-quasiperiodic hyperbolic invariant set F. If, moreover, this set F is attract- 
ing we say that F is a chaotic (strange) attractor. (For definition of hyperbolic 
sets, see Ruelle 1989, Anosov 1995). In this paper, we use only the following 
basic property of hyperbolic sets, so-called Persistence (Ruelle 1989, Anosov 
1995). This means that the hyperbolic sets are, in a sense, stable(robust): if 



(14) generates the hyperbolic set F and 6 is sufficiently small, then dynamics 



(14) also generates another hyperbolic set F. Dynamics (14) and (15) restricted 
to F and F respectively, are topologically orbitally equivalent (on definition of 
this equivalence, see Ruelle 1989, Anosov 1995). It is important to mention that 
a chaos in dissipative systems may be stable, in the sense of structural stabil- 
ity, and although not yet observed in gene networks, structurally stable chaotic 
itineracy is thought to play a functional role in neuroscience (Rabinovitch 1998). 

Therefore, any possible chaotic robust dynamics can be generated by the 
MDF systems, for example, the Smale horseshoes, Anosov flows, the Ruelle- 
Takens-Newhouse chaos, see Newhouse, Ruelle, and Takens, 1971, Smale 1980, 
Ruelle 1989. Some examples of the MDF systems were given in Dancer- Polacik 
1999, Rybakowski 1994, Vakulenko 2000. 

Assertion 2.1 allows us to apply this approach to centralized network dy- 
namics. To this end, assume that (|8| and ^ hold. Moreover, let us assume 

A,; — K Xi, hi = Khi (17) 

where all coefficients hi are uniform in k as k — > 0. We also assume that all 
direct interactions between centers are absent, B = 0. This constraint is not 
essential. 



Since Uj ~ 0{k) for small k, we can use the Taylor expansion for a in (13). 
Then these equations reduce to 

^^ = p,{C,V{v) - h) - Kv, + m{t), (18) 



where pi = ricr'(O), i = 1,2, ...,M and r is a slow rescaling time: r — n^t. Due 
to conditions ( [l7| , the corrections w^ satisfy 

|Wi| < CK. 



Let us focus now our attention to non-perturbed equation ( 18 ) with Wi = 
0. Let us fix the number of centers M. The number of satelhtes N wiU be 
considered as a parameter. 

The next important assertion immediately follows from well known approx- 
imation theorems of the multilayered network theory, see, for example, Barron 
1993, Funahashi and Nakamura 1993. 

Assertion 2.2. Given a number S > 0, an integer M and a vector field 
F = (Fi,...,Fa/) defined on the ball B^' = {\v\ < 1}, F, G C^{B^^), there are 
a number N , an N x M matrix A, an M x N matrix C and coefficients hi, 
where i = 1, 2, ..., N , such that 

\F,{-)-C,W{-)\ciiBM)<5, (19) 

where 

W,{v) ^ a {Mv - h^) , (20) 

where v = {vi, ...,vm) G R*^. 

This assertion gives us a tool to control network dynamics. Assume hi = 0. 



Then equations ( 18 ) with Wi = reduce to the Hopfield-like equations for 



variables Vi = Vi{T) that depend only on r: 

^^KiW{v)^~Xm, (21) 

where I = 1, ..., Af, the matrix K is defined by Kij = piCijRjXj . The param- 



eters V of (21 ) are K, A/, hj and A, 



In this case one can formulate the following result. 

Assertion 2.3. Let us consider a C^-smooth vector field Q{p) defined on a 
ball Bi^ C R^^ and directed strictly inside this ball at the boundary dB^^ : 

F{p)-p<0, pedB'^'. (22) 



Then, for each S > 0, there is a choice of parameters V such that (21) 5 -realizes 



system (14-)- This means that (21) is a MDF system. 



This follows from Assertions 2.1 and 2.2. 



4 A toy model of centralized network 

In this section we consider a simple centralized network that, nonetheless, can 
produce a number of point attractors (stable steady states). Due to its simple 
structure, we can investigate here the robustness of this system. 



Let us consider a central node interacting with many satellites. This motif 
can appear as a subnetwork in a larger scale-free network. In order to study 
robustness, we add noise to the model. We consider two types of stochastic 
perturbations. The first type of perturbations is a Langevin type additive noise 
that can simulate intrinsic stochastic fluctuations of gene expression dynamics. 
The choice of additive noise is for the sake of simplicity, however more general 
multiplicative noise can be used with no change of the results. The second type 
of noise is a shot-like perturbation that can simulate the external contributions 
to noise, caused by the environment. Furthermore, we replace the sigmoid in ([2| 
by a linear function. This is justified in TF - miRNAs networks, where the action 
of satellites (miRNA's) on centers (TF's) is post-transcriptional and produces 
a modulation of the production rate of the center protein. This modulation can 
be modeled by a soft sigmoid or even by a linear function. Moreover, to simplify 
our model, we assume that all satellites are, in a sense, equivalent. 

The network dynamics can be described then by the following equations: 



dui 

dv 
di 
where fi, Q are defined by 



= -Xu, + f,{v)+^,{t), (23) 

= -iyv + Q{u)+^o{t), (24) 



Q{u) = flo + a^Mi, /j(w) = ra{b{v - hi)), 

i=l 

Here ^i are noises, the coefficient A > is a satellite mobility (degradation rate), 
r > is the satellite maximum production rate, b defines a sharpness of center 
action on the satellites, z/ > is a center mobility (degradation rate), a is the 
strength of the satellites feedback action on the center. 

We consider the following type of noises: non-correlated white noise 

{Ut),^,{t'))^p.A,5{t-t') (25) 

where /3i > are intensities, and shot-like noise 

Ut) ^ P^V^5{t - Tj) (26) 

where tj are random shot times following a Poisson process, j3i are noise am- 
plitude coefficients, and rji are random variables distributed uniformly on [0, 1]. 
In numerical simulations we set 5{t — tj) — 1 with a probability po << 1 and 
8(t~- Tj) — with the probability 1 — po, where r, — jSt, St is a time step. Such 
noises ^i can summarize the effect of a strong environment fluctuations on the 
satellite and center expression. 

We study the problem under the following 

Assumption. Let the derivatives of fi and Q satisfy 

fi{v)Q'{u) > for all i,u,v. 



Then one can show, foUowing (Hirsch, 1988) that the dynamics is monotone, 
and, therefore, all trajectories converge to equilibria. The numerical simulations 
confirm this fact. Notice that the above assumption is not needed when satellites 
are fast, because in this case the asymptotic dynamics is one dimensional and 
in dimension one all the attractors are stable steady states (point attractors). 
Although this simple system can not generate chaos or periodic behavior, the 
number of point attractors can be arbitrarily large, and thus this system is 
nonetheless flexible. 

4.1 Multistationarity of the toy model 

Let us fix the signs of the satellite actions on the center assuming that a < 0. 
This restriction is fulfilled in gene networks, where the centers are transcription 
factors (TF) and the satellites are microRNAs (indeed, usually microRNA can 
only repress transcription factors). Let us show that the toy model admits 
coexistence of any number of point attractors. 



Let us make a transformation reducing (231 and (24 1 to a system of two 
equations introducing a new variable Z by 

N N 

Z = ^w„ G(«)=^ra(6(«-/i,)). 



Then, by summarizing eqs. (23), one obtains 



^ = -XZ + Giv), (27) 

dv 

— = -i^v + ao + aZ. (28) 

This system is relatively simple and it can be studied analytically and numeri- 
cally. Since all trajectories are convergent we obtain that the attractor consists 
of equilibria defined by 

XZ + F{Z)=0, F{Z)^G{iy'\aa + aZ)). (29) 

Let P = {b,hi,N,ao,a} be free parameters that can be adjusted. Like to 
the previous section, we can "control" the nonlinearity F hy V and use the fact 
that F{Z) can approximate arbitrary smooth functions. The following assertion 
shows that the system is multi-stationarity with an arbitrary number of point 
attractors: 

Assertion 3.1. Let N be a positive integer. Then there are coefficients 
6, A > 0,r > 0, where i = 1,...,N, v > and hi,ao,a in such a way that 



equation (29] has at least ?i + 1 stable roots that can be placed in any given 
positions in the Z-space. 

The main idea of the proof can be illustrated by Fig. 1 and holds in both 
cases of the Fermi and the Hill sigmoids. Let us make a variable change w = 



\Z/r. The steady states are solutions of the equation rw — F{w), where the 
function F{'w) is close to a step function with N steps; each step is given by the 
function a{'y{w — hi)) that is close to Heaviside step function for large 7. Here 
7 is a parameter that defines the sigmoid sharpness: 

j = abr{i^X)-\ (30) 

The steady states of the system are given by the intersections between the 
graph of F{w) and the straight line of slope r. An elementary argument shows 
that the intersections lying on horizontal segments of the graph of F{w) are 
stable attractors, whereas the intersections on ascending vertical segments cor- 
respond to repellers. 

The position of the i-th step in w-space is hi and its height is r. Under an 
appropriate choice of hi this entails our assertion (see Fig. 1). In the neural 
network theory, 7 is known as gain parameter. This quantity, defined as the 
product of rates on sharpness divided on the product of degradation coefficients, 
gives the maximal possible density of the equilibrium states in w-space. 

It is useful to note that one gets n + 1 attractors on the horizontal segments 
of the step function provided that hi decrease with i. 

Notice that the main condition to obtain flexibility (multistationarity) is the 
sharpness of the sigmoidal function, meaning that the gain parameter 7 should 
be large. The construction is robust: we can vary Wi,b,hi but the number of 
equilibria is conserved. 

4.2 Robustness and stability of attractors 



The roots of Eq.(29) are point attractors and then they are dynamically stable, 
otherwise, they are repellers and unstable. In Fig. 1, attractors correspond 
to intersections of the straight line y — rw with the curve y = F{w), lying 
on horizontal segments of the graph of F. A simple argument suggests that 
the positions of these attractors are robust with respect to variations of the 
thresholds hi . Indeed, a perturbation of hi induces a horizontal shift of the step 
a{j{w — hi)), and the positions of the attractors are only slightly affected. 

More insight into robustness of the centralized toy model can be obtained 
by considering the noisy case £,i ^ 0. 



First, let us consider the case of the Langevin noise (25). We are interested 



in the robustness of the number and positions of the attractors with respect to 



noises ^i{t). Near a point attractor, the equations (27), (28) can be linearized. 



The linearized dynamics is defined by the following matrix H: 

—A /i 
a —V 

where /i = G'{veq). For large 7 and for stable stationary states fj, is small, fj, = 
G'{veq) — ?► as 7 — )• 00. Let us assume that the noises ^i are independent white 
noises. Using standard results from the theory of linear stochastic differential 
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Figure 1: Intersections of the curve y = _F(ui) and the straight line y — rw correspond 
to steady states; intersections on horizontal segments of the graph of F correspond to 
stable steady states. 



equations, see, for example, Keizer 1987, it follows that small deviations 6Z, dv 
from the equilibrium are normally distributed with the density 



p{SZ,Sv) = constexp{-X ■ M^^ • X*''), X = {SZ,Sv) 



(31) 



where M is a symmetric, positively defined, 2x2 covariation matrix with entries 
'7^11, Ti22, '7112 — "7-21 ■ This matrix can be defined by the well known relation 
(the fluctuation-dissipation theorem): 



HM + MH' 



B, 



(32) 



where B — diag{Bf, i3|) and, since the noises are non-correlated, Bi = y X]i=i Pf^ 

B2 = /3o. 

As a result, a characteristic fluctuation amplitude F^ is proport ion al to the 



1/2 .1/2 



} where 6i are eigenvalues of M. Eq. (32) can be 



maximum max{0j^ 

resolved explicitly and Oi can be found. 

Now we can investigate the following problem: how to tune the parameters 
A, V and a to obtain the minimal fluctuation amplitude Fa with respect to the 
noise under a given multistationarity level (this means 7 = 70 >> 1 is fixed 
but we can vary the degradation rates v,X). This optimization problem can 
be resolved numerically. The results, which describe the optimal Xopt,i^opt as 
functions of 7, are as follows. 

The case A), B2 >> i?i, the center is under a stronger noise than the 
satellites. Then the center degradation rate Vopt should be large, and Xopt is a 
small, decreasing in 7 function. 
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The case B), i32 < -Bi, the center is under smaller noise than the satellites. 
Then the center degradation rate Vopt should be smaller, Xopt > Vopt, and the 
both parameters are decreasing in 7. This situation is illustrated by Fig. 2. 
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Figure 2: Optimal degradation parameters Vopt, ^opt (minimizing the eigenvalues of 
the matrix H) as functions of the gain parameter 7 in the case B2 < B\, when the 
center is under smaller noise than the satellites. 

The classical ideas of the invariant manifold theory, discussed in the preced- 
ing section, allow us to systematize these results. The centralized network can 
function under two main and quite opposite regimes. The first one arises when 
X» V. Then the satellite dynamics is slaved by the center motion. The center 
dominates and such a regime can be named power of the center. This regime 
is stable if /Sq is small, but Pi are large (the noises act on satellites mainly, 
case B). Considering that the noise intensity is larger for those components 
that are expressed in larger copy numbers, the case should be representative 
for miRNA-TF networks, when miRNA are in smaller copy numbers than the 
transcription factors. In this case the noise perturb satellite states {ui) but, 
since the satellites are controlled by the center state w, satellites return to the 
normal states and dynamics is robust, the noise does not damage the attractor. 
The opposite regime is when X « v. Then, opposite to the previous situation, 
the center dynamics is slaved by the satellites motion. Such a regime can be 
named satellite democracy. This regime is stable when /3o is large, but j3i are 
small (the noise acts stronger on the center, case A). Here the noise can perturb 
the center state (w) but this state can be restored by satellites. The large time 
dynamics is robust, again the noise does not damage the attractor. 

Similar results are illustrated in the case of a shot noise in Fig. 3. 

So, we obtain an interesting connection between robustness, multistation- 
arity and network rate: to support robustness and multistationarity in a noisy 
situation, we should decrease the degradation constants. Multistationarity of 
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molecular switches is important in decision making processes in differentiation, 
development, and immune response of the organisms. Our finding means that 
noise protected switches are necessarily slow. 
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Figure 3: Numerical simulations of the system's trajectories under shot noise. The 
parameters were as follows: N — 6, (j{z) — (th(z) with p — 4, b = 20, Ka = 1, 
/li = i, r = 1, A = 5, ao = 0.3u, x G [0, 10] and t G [0, 120]. The parameter A = 
for all i beside i = 3, where Ps = 50. The two functioning regimes correspond to 
different values of z/, a, namely, a — 50, i/ = 5 (satellites democracy (SD), a)), and 
a = 5, 1/ — 0.5 (power of the center (PC), b)). In the both cases the system shows 
multistationarity. For the chosen initial data trajectories converge to an attractor 
V ~ 5.3 as t >> 1 in the PC regime and also in the SD regime v — >■ 3.3. This means 
that the fast center loses the attractor control, while the slow center controls dynamics 
even under large deviations. 



5 Conclusion 

We have considered networks with two types of nodes. The w-nodes, called 
centers, are hyperconnected and interact one to another via many u-nodes, 
called satellites. We show, by recently advanced mathematical methods, that 
this centralized network architecture, allows us to control network dynamics 
to create complicated dynamical regimes. This network organization creates 
feedback loops that are capable to generate practically all kinds of dynamics, 
chaotic or periodic, or having a number of equilibrium states. This strong 
flexibility could be crucial for adaptive biological functions of these networks. 

Using the simple example of a motif with a single center, we also argued that 
centralized networks can perform trade-offs between flexibility and robustness. 
To support both flexibility and robustness in a noisy situation, the network 
should function in a slow manner, i.e, we propose slow-down as a way to increase 
stability. 
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Which of the nodes should be slowed-down depends on the fluctuations. 
Basic ideas from the invariant manifold theory show that if the noises act on the 
satellites, then, in order to conserve dynamics and the attractor structure, the 
center should be slow and controls the satellites (we called this regime power 
of the center). In the opposite case, when the noise acts on the center, the 
satellites should be slow in order to control the center and the global dynamics 
(we called this regime satellites democracy). 

We did not consider here extrinsic noise or parametric variability of the 
system, that we plan to study in the future. We also think that the slow-down 
effect could be observed in all systems where there is a separation into slow and 
fast variables, independently of architecture. 
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